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Abstract. It is shown that the polynomial 

p{t) = Tr[{A + tS)"] 

has positive coefficients when m = 6 and A and B are any two 3-by-3 complex 
Hermitian positive definite matrices. This case is the first that is not cov- 
ered by prior, general results. This problem arises from a conjecture raised by 
Bessis, Moussa and Villani in connection with a long-standing problem in the- 
oretical physics. The full conjecture, as shown recently by Lieb and Seiringer, 
is equivalent to p{t) having positive coefficients for any m and any two n-by-n 
positive definite matrices. We show that, generally, the question in the real 
case reduces to that of singular A and B, and this is a key part of our proof. 



1. Introduction 

In P^, while studying partition functions of quantum mechanical systems, a 
conjecture was made regarding a positivity property of traces of matrices. If this 
property holds, explicit error bounds in a sequence of Fade approximants follow. 
Recently, in [5], and as previously communicated to us [1], the conjecture of [T] was 
reformulated as a question about the traces of certain sums of words in two positive 
definite matrices. 

Conjecture 1.1 (BMV). The polynomial p{t) = Tr [(A + tS)™] has all positive 
coefficients whenever A and B are n-by-n positive definite (PD) matrices. 

The coefficient of t^ in p{t) is the trace of Sm,k{A, B), the sum of all words of 
length m in A and B, in which fc S's appear (sometimes called the fc-th Hurwitz 
product of A and B). In [3], among other things, it was noted that, for m < 6, 
each constituent word in Sm,kiA, B) has positive trace. Thus, the above conjecture 
is valid for m < 6 and arbitrary positive integers n. It was also noted in [J that 
the conjecture is valid for arbitrary m and n < 3. Thus, the first case in which 
prior methods do not apply and the conjecture is in doubt, is m = 6 and n — i. 
Even in this case, all coefficients, except Tr[5'6,3(A, i?)], are known to be positive 
(also as shown in [3]). Our purpose here is to show that the remaining coefficient 
Tr[5'6,3(A, J5)] is nonnegative when A and B are 3-by-3 positive definite matrices, 
which requires notably different methods (some summands of S'6.3(A, _B) can have 
negative trace [4]). It follows that the conjecture is valid for m = 6, n = 3, our new 
result. A key tool is that it suffices to prove the conjecture for singular (positive 
semidefinite) matrices. 

The coefficients Sm,k{A, B) may be generated via the recurrence: 

{A,B)^Sra,k{A,B)B + {A,B)A 
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(variants are available). The following lemma will be useful for computing the Sm.k- 
We give an algebraic proof although a purely combinatorial proof is also available. 

Lemma 1.2. For any two n-by-n matrices A and B, we have 
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Since Sm,k{Ay, B) = y™ ^ Sm.k{A, B), it follows that the coefficient oi t^ in the 
last expression above is just 

mTr[v45„_i,fc(A,S)] - (m - A:)Tr[5„^fc(A, B)], 

which proves the lemma. □ 



2. Reduction to the Singular Case 

Of course, when A and B are Hermitian, Sm,kiA, B) is Hermitian, but even 
when A and B are n-by-n real symmetric PD matrices, n > 2, Sm,kiA, B) need 
not be PD. Examples are easily generated, and computational experiments suggest 
that it is usually not PD. We want to show that Tr[SQ^3{A, B)] is nonnegative for 
3-by-3 positive definite A, B. This is subtle as 5*6,3 (A, i?) need not have positive 
eigenvalues, and as some words within the 5*6,3 (A, B) expression can have negative 
trace [41. A main component of our argument is based on the following technical 
observation. 

Theorem 2.1. Let B be any real n-by-n matrix, and let A — diai,g{l,xi, ... ,Xn-i). 
Suppose that a — (ai, . . . ,a„_i) £ [0, 1]"^^, and let D = diag(l,c?i, . . . be 
such that di = if Oi — 0, and di = I otherwise. If a achieves the minimum of the 
function f : [0,1]"""'^ — » R given by f{xi,...,Xn-i) = Tr[Sm,k{A, B)], then, with 
A' — diag(l, fli, . . . , a„_i), we have 



/(ai, ■ 



_i) =Tr[5„,,fe(A',B)] - 



m — k 

Proof. Let A', B, D, and a = (ai, . . . , a„_i) G [0, be as in the hypotheses 
of the theorem. First suppose that A' = D. Then, it is clear that the formula in 
the theorem reduces to the identity in Lemma 11.21 When A' ^ D, consider the 
differentiable function g : [—1/2, 1] ^ R given by 

g{z) = Tr SraM ("f^j""'^ 
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By hypothesis, a G [0, 1]" ^ achieves the minimum for /. Consequently, it follows 
(from basic variational techniques) that 

(2.1) 



Next, notice that, 
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In particular, at z = 0, the above expression evaluates to 



Tr 



J2 + tBY'^ {D - A') {A' + tBf 



= Tr 
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Finally, observe that Sm,k{A'y,B) ~ y™ Sm.k{A' , B) so that the coefficient of 



t" in ([231) is 

mTi[DS,n-iMA', B)] - (m - k)Tr[S^,k{A' , B)]. 
It follows, therefore, from (|2.1[) that 

777 

Tr[5™,fc(^', B)] -Tr [D5„,_i,fe(A', B)] . 

m — K 

This completes the proof. □ 

Example 2.2. As an example of the theorem, let m — A, n — 3, k — 2, and 

"10 

, A= xi 

_ X2 

A straightforward computation gives us that 

Ti[S4,2{A, B)] = 20-4x1 +8 a;!^ - 12xiX2 +42x2^ 

Tr[53,2(A5)] = 9 + 18x2. 
The minimum of Tt[S4,2{A, B)] is achieved by xi = 7/25, X2 = 1/25, and one has 

Tr[54,2(A',S)] = 2Tr[53,2(A',i3)] = ^. 

25 

Let A, B, and / be as in Theorem l2.1l If we are fortunate enough that / achieves 
a minimum /(a) with a G (0, 1]"~^, then D is the identity matrix and the theorem 
statement simplifies to the following. 
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Corollary 2.3. Suppose that f as in Theorem \2.1\ achieves a minimum /(a) with 
a G (0,1]"^"'^. Then, the nonnegativity of TT[Sm-i,kiA' , B)] implies the nonnega- 
tivity ofTT[SmAA',B)]. 

To see the importance of this corollary, we next examine the real version of 
Coniecture ll.il Suppose we know that the conjecture is true for the power m — 1 
and also suppose (by way of contradiction) that there exist n-by-n real positive 
definite matrices A and B such that Tr[S',„,fe(^, B)] is negative. Then, in particular, 
(by homogeneity) there are real positive definite A and B with norm 1 such that 
Tr[S'm,A:(^, B)] is negative (here, we use the spectral norm [6l p. 295] so that 
for positive semidefinite A, it is just the largest eigenvalue of A). Let M be the 
(compact) set of real positive semidefinite matrices with norm 1 and choose {A, B) S 
M X M that minimizes Tr[S'„i_fc(A, _B)]; our goal is to show that this minimmn is 0. 
By a uniform (real) unitary similarity we may assume that A — diag(l,ai, . . . ,a„_i) 
is diagonal with 1 > ai > • • • > a„_i > 0. 

Corollary 12.31 then tells us that A must be singular, because by induction, 
Ti[Sm-i,k{A^B)] will be nonnegative for all positive semidefinite A and B. By 
symmetry, it also follows that B is singular. We combine these observations into 
the following theorem. 

Theorem 2.4. Suppose that Tr \{A + ii?)™^^] has all positive coefficients for each 
pair of n-by-n real positive definite matrices A and B. Ifp(t) — Tr [{A + tB)"^] has 
all positive coefficients whenever A, B ^ are singular n-by-n real positive definite 
matrices, then p{t) has all positive coefficients whenever A and B are arbitrary 
n-by-n real positive definite matrices. 

3. Symbolic Real Algebraic Geometry 

In this section, we discuss the symbolic algebra preliminaries necessary for solving 
the m = 6, n = 3 case of Conjecture 11.11 Let R = Q[a;i, . . . , a;„], and let /, J be 
two ideals of R. The quotient ideal of / by J is the ideal of R given by [21 p. 23] 

{I : J) ^ {f e R: fg e I for all g £ J}. 

We can iterate this process to get the increasing sequence of ideals 

/ C (/ : J) C (/ : J2) C (/ : J3) C • • ■ . 

This sequence stabilizes to an ideal called the saturation of / with respect to J (see 
B P- 15]): 

(/, J°°) = {f eR:3meN with /™ • J C /}. 

If / is any ideal in R, let ¥{1} denote the set, 

Vil) = {(ai, . . . , a„) e C" : /(ai, . . . , a„) = for all / e /}. 

From these definitions, it is easily verified that for any two ideals, /, J C i?, 

ViI)\ViJ) C Vil : J°°). 

For our particular application, we will be interested in proving that V{I) \ V{ J) 
contains no elements in (0, 1)". Let P denote the saturation ideal (/ : J°°). If we 
are fortunate enough to find that P = (1) = Q[a;i, . . . , x„], then there are no points 
in V{I) \ V{J) (and hence none in (0, 1)"). One difficulty with this approach is that 
these new saturations do not always produce unit ideals. One more idea is needed, 
which we describe below. 
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If K is an ideal of i?, the elimination ideal 2, p. 25] of K with respect to Xi is 
Ki = Q[a;i]ni?. The Xi-coordinates of elements in V{K) are elements in V{Ki). For 
our purposes, we need only verify that for a saturation P, there is an elimination 
ideal Pi of P such that V{Pi) contains no numbers in (0, 1). 

Normally, a procedure such as the one outlined above would be relatively in- 
tractable (the symbolic algorithms are doubly exponential in nature). Our reduc- 
tions give us enough efficiency to complete a proof computationally. We performed 
our computations using the symbolic algebra system Macaulay 2. 

4. The Case m = 6, n = 3 

The remainder of this article is devoted to a technical consideration of the case 
rn = 6, fc = 3, n = 3 which is the content of the theorem below. 

Theorem 4.1. The polynomial p{t) ~ Tr[(yl + i_B)'™] has positive coefficients when 
m = 6 and A and B are any two 3-by-3 positive definite matrices. 

Proof. Suppose that there exist 3-by-3 (complex Hermitian) positive definite ma- 
trices A and B such that Tr[iS'6,3(A, B)] is negative; we will derive a contradiction. 
Performing a uniform unitary similarity and using homogeneity, we may assume 
that A and B are of the form, 
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in which 1 > r > s, a, 6, c > 0, and x,y, z G C. If x,y, z > 0, then we clearly have 
a contradiction. Otherwise, perform a simultaneous diagonal unitary similarity on 
A and B (a similarity by a diagonal matrix with entries on the unit disc) making 
x,y > 0. This does not change the trace of Se,3{A, B). 

We next show that we may assume z G M. A computation of Tt:[Sq^3,{A,B)] 
reveals that it has the form w = azl + /3z + ^z + 6, in which a, /3,j,S > 0. Since 
w is real, we have 

w — Re{w) — azz + /3Re(z) + 7Rc(2) + S 

> aRe{z f + (3Re{z) + jRe{z) + 5. 

Consequently, it follows that we can assume z is real and negative. Theorem 12.41 
now applies, so that it is enough to verify the claim with s = and det{B) = 0. 

Since B is positive semidefinite, we have a6 — a;^ > 0. If 6 = 0, then a; = 0, and 
an easy computation shows that 

Tr[S'6,3 {A, B)] =6z'^c + 24 az^ + 20 + 6 r^y'^c > 0, 

a contradiction. Therefore, we must have 6 > 0. A similar computation also shows 
that a,x,y > 0. 

Next, we prove that c > 0. Since 

det(i?) = 2 xzy + abc — ay^ — x^c — z'^b — 0, 

it follows that when c = 0, we have 2a;?;z = bz'^ + ay^. From this, it is clear that 
z < is impossible, and therefore z = 0, a contradiction. Finally, if ab = x^, then 
from det(i3) = 0, we have that 2xyz = 6z^ + x'^y'^/b. This implies again that z — 0, 
another impossibility. Hence, ab — x^ > 0. 
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Summarizing these observations, we may assume that 

X —z 



B 



b 

X b y 

xbzy+z^b'^ 



(4.2) D = 



in which u,b,x,y > and z > 0. Furthermore, if r = 1 or r = 0, then [SJ Theorem 
4] (along with a straightforward continuity argument) imphes that TtISq^s^A, B)] 
is nonncgative. Therefore, we may assume that < r < 1. 

A direct computation shows that 6^M^Tr[S'6,3(A, B)] is a polynomial p(r, x, y, z, u, b) 
€ Z[r, X, y, z, u, b]. The negative terms in p factor as 

(4.1) - 12 b^u^xzy {r^ + r+l) . 

We shall verify that the minimum of p{r,x,y, z,u,b) over r,x,y, z,u,b G [0,1] is 
0, which will prove the claim (by homogeneity of the matrix B in the variables 
x,y,z,u,b). 

If any of x, y, z, u, or b is zero, then we are done by (j4.ip : therefore, we begin by 
determining the critical points of p in (0, oo)^. This amounts to a calculation of 

'dp dp dp dp dp dp^ 
^dr' dx' dy' dz' du db ^ 

which is an ideal in the ring Q[r, x, y, z, u, h\. We are interested in verifying that the 
set of points V{D) \ V{rxyzub) contains no element in (0, 1)^. From the discussion 
in the previous section, it suffices to verify this claim for V{D : {rxyzub)°°). 

Let P — {D : {rxyzub)°°). Using Macaulay 2, it can be checked that P is the 
unit ideal Q[r,x,y,z,u,b\. It follows that the minimum of the function p above 
must occur when one of the x, y, z, u, 6 is 1 (in other words, on the "boundary"). 

This process now continues, recursively, by next finding the critical points of 
the functions p{r, 1, y, z, m, b), . . . ,p{r, x, y, z, u, 1), and checking that they either do 
not occur in (0, 1)''' or that the function is nonnegative when they do. As noted 
before, a difficulty is that these new saturations do not always produce unit ideals. 
Therefore, we finish by showing that for each saturation P, there is an elimination 
ideal Pi of P such that V{Pi) contains no positive numbers in (0, 1). Since each 
Pi is generated by a single-variable polynomial, we use Sturm's algorithm to verify 
such a claim symbolically. These computations were also performed in Macaulay 
2. This completes the proof of the theorem. □ 

As a final remark, we should note that there are some good tools for the nu- 
merical exploration of such problems. Namely, the program SOSTOOLS written 
by Prajna, Papachristodoulou, and Parrilo is an excellent resource for investigating 
real algebraic systems Q 
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